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ABSTRACT 

We present Spitzer observations of the unusual variable V838 Monocerotis. Extended emission is 
detected around the object at 24, 70 and 160 /im. The extended infrared emission is strongly correlated 
spatially with the HST optical light echo images taken at a similar epoch. We attribute this diffuse 
nebulosity to be from an infrared light echo caused by reprocessed thermal emission from dust heated 
by the outward-propagating radiation from the 2002 eruption. The detection of an IR light echo 
provides an opportunity to estimate the mass in dust of the echo material and hence constrain its 
origin. We estimate the dust mass of the light echo to be on the order of a solar mass - thereby 
implying the total gas plus dust mass to be considerably more - too massive for the echo material to 
be the ejecta from previous out burst/mass- losing events. This is therefore suggestive that a significant 
fraction of the matter seen through the light echo is interstellar in origin. Unresolved emission at 24 
and 70 /im is also seen at the position of the central star possibly indicating the presence of hot dust 
freshly condensed in the outburst ejecta. 

Subject headings: infrared: stars-novae, cataclysmic variables - stars: individual (V838 Monocerotis) 



1. INTRODUCTION 

The eruption of V838 Monocerotis (V838 Mon) in Jan- 
uary 6, 2002 has introduced a seemingly new kind of ob- 
ject in the real m of cataclysm ic variables. It was detected 
in eruption by iBrown 1 l)2002j) with a peak outburst am- 
plitude of Vmax = 10 which was followed by two more out- 
bursts - within the next two months - reaching 6.7 and 
7 magnitudes respectively. V838 Mon also showed fast 
cooling, on a timescale of a few months, to a cool, late M 
spectral type or beyond ijEvans et al.ll2003T) . The multi- 
peaked outbursts in the object and the decrease of its ef- 
fective temperature with time suggested that the object 
was different from a classical nova or other known classes 
of eruptive variables (Munari et al. 2002; Kimeswenger 
et al. 2002). A n expanding light-e cho was also seen 
around the star l)Henden et alJ 12002') whose expansion 
and accompanying m orphological chang es are strikingly 
illust rated bv HST llBpnd et alJ 127)031) and other im- 
ages l|Crause et al.ll2005|) . Distance estimates to the ob- 
ject, converging to a large value in the range of 7-10 
kpc, have been made based on varied approac hes such 
as the rate of the expa nding light echo (e.g., iTvlendal 
20041 iCrause et al.l2005|) . a nd the detection of Si O maser 
emission from the source l|Desfuchi et all 120051) . Con- 
siderable near-infrared studies of the source have been 
done yielding spectra showing an oxygen rich atmosphere 
with prominent molecular features of CO, AlO , SiO, 
TiO and water llLvnch et all 120041 lEyans et alJ l2003t 
iBaneriee fc Ashokll2002t iRushton et al.ll2005|) ~hough 
not firmly established, similarities in outburst properties 
between V838 Mon and two other potential analogs viz. 
V4332 Sgr and M31-RV, suggest that they could be uni- 
fied into a new class of eruptive variables. While the 
cause for the intriguing outbursts in such objects is yet 
to be securely established, various mechanisms have been 
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proposed viz. the merg er between main sequence stars 
( Soker & Tvlcnda 2003 ); the capture of multiple plane- 
tary companions ijRetter & : Marom|[2003') or a late flash 
in a born-again ACB star l)Lawlorll20 05'). Aspects relat- 
ing to the spectral type of the progenitor and whether it 
is a single or binary star have also been the subject o f 
investigations ijMunari et al.ll2005t iTvlenda et alJ l2005). 

Here, we present our results on V838 Mon from GO 
Cycle 1 observations of the Spitzer Space Telescope 
(Spitzer). The highlight of the Spitzer observations is the 
striking infrared echo seen around V838 Mon in the mid- 
and far-IR images. While optical light echoes around 
novae or supernovae are rare but not unknown, an in- 
frared "echo" is a rare complementary phenomenon - 
we are only aware of on e other instance of a resolved in- 
frared light echo, Cas A ijKrause et al.l2005T) . When pho- 
tons from the illuminating source interact with the dust 
grains of the echo material, they can either be absorbed 
or scattered according to the albedo of the grain. Pho- 
tons scattered into the line of sight result in an optical 
light echo, a direct "image" of the impinging radiation. 
On the other hand, IR "echoes" result from the absorp- 
tion of the impinging radiation; the thermalized energy 
of the absorbed photons is re-emitted in the IR. Our ob- 
servations are presented in SJ2I Details of the nature of 
the IR emission are discussed in ^3 In 21 the important 
question regarding the origin of the light echo material 
around V838 Mon is discussed. 

2. OBSERVATIONS AND DATA REDUCTION 

V838 Mon was imaged with the Multiband I maging 
Photometer for Spitzer (MIPS; iRieke et al]l27)04j) at 24, 
70 and 160 /im in 2004 and 2005. The log of the obser- 
vations is given in Tabled 

The MIPS d ata were reduced us ing the Data Analy- 
sis Tool (DAT; iGordon et all 12005]) . Images at 24, 70, 
and 160 /im are shown in the first column of Figure 1. 
While at 24 /im, V838 Mon appears as a bright point 
source (hard saturation in the core), extended nebulos- 
ity offset from the central star is readily apparent at 70 
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TABLE 1 

Log of Spitzer Observations for V838 Mon 





Instrument 


A 


AA 




Exp. per frame 




Integration 


AOR 


Date 


Array 


/tin 


jum 


Mode 


(sec) 


Cycles 


on-source (sec) 


Key 


2004 Oct 14 


MIPS/24 


23.7 


4.7 


Fixed Photometry 


3 


2 


92.3 


10522624 


2005 Apr 01 


MIPS/70F 


71 


19 


Phot, w/ 8 clus. pointings 


10 


2 


615.3 


10523648 


2004 Oct 15 


MIPS/160 


156 


35 


Phot, w/ 4 clus. pointings 


3 


2 


50.2 


10523904 



TABLE 2 
Observed Flux Densities 



X eff 




Flux Density (Jy) 1 






Total a 


Core b 


Extended 


23.68 


15.97 


15.06 


0.91 


71.42 


14.72 


3.82 


10.9 


155.89 


17.53 




17.53 



Absolute calibration errors are 5, 10, and 10%, respectively. 
a Total flux density, 80" aperture. 
b Unresolved point source from PSF fitting. 
c Residual emission (Total— PSF). 



and 160 /im. Detailed comparison of the 24 /im radial 
profiles between the target and an observed standard star 
indicates that the profiles matched very well between the 
first bright and dark Airy rings, suggesting most of the 
flux at the central region is from an unresolved source, 
and the extended nebulosity contributes very little flux 
at the core. To correct the saturation in the 24 /im im- 
age, an observed point spread function (PSF) was scaled 
to match the brightness of the first bright Airy ring. The 
saturated core of V838 Mon was then replaced with the 
scaled PSF. Total fluxes observed through 80" apertures 
are reported in column 2 of Tabled To separate the ex- 
tended emission from the unresolved point source, PSF 
subtraction was performed on the 24 and 70 /im images 
(at 160 /im, there is no discernible unresolved central 
source) . The flux density of the central source estimated 
from the subtracted PSF is reported in column 3 of Table 
[21 The PSF subtracted images are presented in column 
2 of Figure 1. Extended emission is clearly evident at 
24 /im in the PSF subtracted images, especially a bright 
extension to the south of V838 Mon. The total diffuse 
flux observed toward V838 Mon, obtained by subtracting 
the unresolved point source flux from that measured in 
the 80" aperture, is given in column 4 of Table [3 

In order to compare our infrared images with the opti- 
cal light echo, we have also obtained archival HST ACS 
images of V838 Mon (Hubble Heritage, GO/DD 10392). 
The data were obtained with the ACS/WFC using the 
F435W, F606W, and F814W filters and were completed 
on 23 Oct 2004, within days of our 24 and 160 /im ob- 
servations and a few months of our 70 /im observations. 
The F814W image convolved with the MIPS beams is 
shown in column 3 of Figure 1 for comparison. A three 
color composite image was constructed using all three 
HST filters and is presented in Figure 2 with our 24, 70, 
and 160 /im contours superposed. The 24—160 /im spec- 
tral energy distribution (SED) of the extended nebulosity 



around V838 Mon is shown in Figure 3. 

3. SOURCE OF THE IR EMISSION 

The infrared emission associated with V838 Mon con- 
sists of an extended component associated with the light 
echo and a compact, unresolved component spatially co- 
incident with the central star. The unresolved compo- 
nent evident at 24 and 70 /im possibly indicates the con- 
densation of newly formed dust around the central star of 
V838 Mon. Evidence that such a dust shell had formed 
after th e outburst has alrea dy been seen from the mid-IR 
data of ILvnch et all (|2004|) in early 2003. These authors 
show the prominent presence of a 650K component in 
their SED in the 8- 13 /im region. The det ection of water 
around V838 Mon ijBaneriee et all 12005). from near-IR 
data in 2002-2003, also supports the existence of a cool 
envelope at ~ 800K. Our MIPS data for the unresolved 
component (Tabled Col. 3) are consistent with a tem- 
perature T>100 K. However, a more detailed and accu- 
rate characterization of the properties of this unresolved 
component - such as determining its temperature more 
accurately, estimating the mass of the newly formed dust 
around the central star to get a better insight on how 
much mass is lost in V838 Mon-type of outbursts, and 
also a detailed study of its spectrum in the 5-40 /im re- 
gion - will be undertaken in a separate study involving 
additional Spitzer data. 

The extended component of emission in V838 Mon 
is strongly indicative of an "infrared light echo". The 
strong spatial correlation of the infrared emission with 
the contemporaneous optical light echo (Oct 2004; see 
Figs. 1 & 2) indicates that the IR emission is also linked 
to the radiation from the initial outburst. The detection 
of light echo emission at wavelengths where the scattering 
efficiency of typical dust grains is small relative to their 
absorption efficiency indicates that the IR light echo - 
unlike the optical light echo - is very unlikely to be due 
to scattering, but rather thermal emission from grains 
heated by the outburst pulse. 

Simple energy arguments also support a thermal ori- 
gin for the IR emission. During the first two months 
of its multi-peaked outburst, when V838 Mon can be 
thought to have emitted most of its radiation in the form 
of a prolonged pulse, the star showed an SED well rep- 
res ented by effective te mperatures between 5250 to 4250 
K l|Munari et al.ll2002l) . The V magnitudes ranged from 
10 to 6 during the outburst, corresponding to an en- 
ergy flux output in the V band, corrected for redden- 
ing using E(B - V) = 0.5, of between 3.9xl0" 16 and 
1.6xl0 -14 W cm~ 2 /im -1 . Assuming a black body tem- 
perature of 5000 K and extrapolating to 160 /im, this 
results in ^10~ 24 to 10~ 23 W cm~ 2 /im -1 available from 
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Fig. 3.— The 24-160 /urn SED of the extended component (light 
echo) of V838 Mon in late 2004 and early 2005. The filled circles are 
the observed flux densities of the light-echo data (Table 2, column 
4) which were fit by two modified blackbodies with temperatures 
of ~25 and ~63K (dashed and dotted lines respectively). The 
resultant sum of these fits is shown by the continuous line. See the 
text for further details. 



the outburst in the 160 /im band. However, our ob- 
served flux density at 160 /mi of 17.53 Jy corresponds 
to an energy flux of ^2.0xl0 -19 W cm -2 /im -1 . Even if 
we unrealistically assume 100% scattering efficiency, the 
outburst provides 3 to 5 orders of magnitude too little 
energy to power the 160 /im emission through scattering 
alone. Thus scattering is ruled out as the source of the 
infrared light echo. Instead, the heating of dust grains by 
optical and UV photons from the outburst pulse and the 
subsequent re-processing of the absorbed energy into the 
thermal infrared remains the likely source of the infrared 
light echo in V838 Mon. 

4. ORIGIN OF THE IR ECHO MATERIAL: INTERSTELLAR 
OR CIRCUMSTELLAR? 

The origin of the light echo material h as been the 
source of so me debate in t h e liter ature. iBond et aD 
(2003) and Ivan Loon et alJ l)2004|) suggest that the 
echo material has been produced in previous mass loss 
episodes while, eg. iTvlendal l)2004D has argued that the 
light echo material is interstellar in origin. This is an 
important aspect to address since it can help estab- 
lish whether the progenitor has undergone previous out - 
bursts (eg., AGB-like behavior; Ivan Loon et al.l {2004)) 
or steady mass- loss through winds in the past, and there- 
fore give a deeper insight into the nature of the central 
star. The discovery of an IR light echo provides an oppor- 
tunity to estimate the mass in dust of the echo material 
and hence constrain its origin. The flux density mea- 
sured for the light echo material is given in the fourth 
column of Tabic 

The infrared flux received from a mass Md of dust 
grains can be estimated from 

Frfust (A) = ^^K i (A)S A (T A ,) (1) 

i 

where D is the object's distance, k\ is the wavelength 
dependent mass absorption coefficient, BxiTp) is the 
Planck function at dust temperature Tq, and the sum ex- 
tends over all dust components contributing to the emis- 
sion at the observed wavelengths. To estimate the mass 



of the emitting echo material observed by MIPS, we can 
use Eq. ^and our measured MIPS flux densities. While 
the echo likely arises from material with a range of tem- 
peratures, a minimum of two components is required to 
fit the 24-160 /im echo emission and we therefore elect 
to fit the data with this minimum number of compo- 
nents consistent with the data. We use Eq. combined 
with a Monte Carlo technique to explore the range of 
temperatures and masses that can fit the data. Allow- 
ing a range of size from 0.1 — 1 /im an d compositions 
(silica te /a morphous carbon; i ndices from lLaor &: Draind 
(1993) and lZubko et al.l<|1996l) . respectively), the derived 
temperatures and masses for the two components are 
~ 63lf °/ 25 -5 K and °- 6 -o.55 x 10 -3 /1.6l^ M Q , respec- 
tively, assuming a distance of 8 kpc. Estimates of the un- 
certainty are derived from the % 2 distribution. The tem- 
perature and mass uncertainty in the warm component is 
due largely to the lack of constraining short wavelength 
data. However, it is important to note that the domi- 
nant component of the derived mass of the light echo, 
arises from the cold 25K component which is reasonably 
well constrained in the fit in Figure 3. Even assuming a 
lower limit on the distance to V838 Mon of 5 kpc, the 
mass in the cold component must be > 0.2 M a . As can 
be seen in Figure 3, the need for a large mass of colder 
dust is primarily driven by the 160 /tm point. The close 
correspondence of the distribution of emission at 160 /im 
and the optical light echo (see the bottom row of Fig- 
ure 1 and the last panel of Figure 2) argue that the cold 
emission is associated with the echo material; indeed, the 
emission at 160 /im rapidly falls off away from the echo 
location. Thus the mass we are estimating from our IR 
flux measurements is directly associated with the echo 
material. With a gas-to-dust ratio of 100, we estimate 
the total mass of the material in the echo to be on order 
of a few tens to a few hundred solar masses. At a dis- 
tance of 8 kpc, the size of the IR echo is ^3 pc; a large 
diffuse interstellar cloud of this size may have a mass of 
a few tens of M Q while dens er clouds could range up to 
a few hundred M in mass (|Whittetl 12003]) . Therefore, 
the observed mass of the echo material is consistent with 
an origin in interstellar material. 

Conversely, it is extremely unlikely that a single pre- 
vious outburst could have ejected a total mass of the or- 
der of 10 M or above. In novae eruptions, typically 
10 -5 Mq matter is ejected; an amount in the range 
10~ 2 — 10~ 5 Mq is suggested from limits /estimates of 
the ejecta mass in V838 Mon (|Baneriee fc Ashokll200l 
iRushton et al.l 12003^1 . Steady mass loss from the pro- 
genitor in the past also does not appear to be favored 
as the sole contributor to the echo matter. The radius 
of the emitting zone (i?), either derived from the light 
travel distance since the outburst ~ 2.75 years ago or 
from the angular radius of the echo in December 2004 
(70"), is on the order of parsecs. The amount of mat- 
ter that could have been lost in such a zone of radius 
R from a star with a wind velocity v, ejecting mass at 
a rate of 10 _n M /year can then be easily computed. 
Even a massiv e AGB star with a hi gh mass-loss rate of 
10~ 5 M / year l)van Loon et al.120 05) in a slow wind with 
v = 10 km/s can only generate a total echo mass (gas 
plus dust) of ~1 Mq. Even this is just at the margin 
of the lower limits of our estimate for the echo material 
mass and we are thus led to believe that a considerable 
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part of the echo material is likely to lie in an intervening 
dust sheet. This is consistent with some of the models 
explaining the light echo expansion rat e and observe d 
morphological changes in the light echo l)Tvlendal l2004'). 
However our calculations do not completely preclude a 
component of circumstellar origin for the echo material; 
indeed the presence of an unresolved core of emission at 
24 and 70 /im argues that some of the material may well 
be "circumstellar" in origin. 

5. CONCLUSION 

We report the detection of a rare infrared "light echo" 
around the unusual variable V838 Mon. Our Spitzer 
data reveal both an unresolved hot component and an 
extended cooler component. The unresolved component 
is attributed to the formation of dust in the ejecta of the 
2002 outburst. The extended component is spatially co- 
incident with the optical light echo and is likely caused by 



thermal emission from dust heated by the energy of the 
outburst pulse. The large derived mass of the infrared 
light echo makes it unlikely that the echo material is the 
remnant of previous mass loss episodes in V838 Mon, 
but rather is interstellar material along the line of sight 
to V838 Mon. 
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Fig. 1. — Multi-wavelength V838 Mon images. Column 1: MIPS 24, 70, and 160 /an images (top to bottom); Column 2: PSF subtracted 
images (no point source detected at 160 /im, see text); Column 3: HST F814W image convolved to the MIPS resolution at 24, 70, and 
160 (im images (top to bottom) after removal of field stars in the HST image. MIPS beam sizes are indicated by the white circles with 
FWHM of 6", 18", and 40". All images are displayed in a field of view of 160" by 160", with N up and E toward the left, and in false-color 
logarithmic scale with brightness and contrast adjusted for best presentation of each image. The displayed surface brightness (top to 
bottom) ranges from 10 — 2 to 10 2 mjy/arcsec 2 at 24 /im, 5xl0 — 1 to 8 mjy/arcsec 2 at 70 /im, 6xl0 _1 to 2 mjy/arcsec 2 at 160 /im for the 
first column; and from 3xl0 — 2 to 5xl0 — 1 mjy/arcsec 2 at 24 /im, and 3xl0 — 1 to 4 mjy/arcsec 2 at 70 /im for the second column. 




Fig. 2. — Morphological comparison between the optical and infrared light echo. The background is the three color composite HST image from 23 Oct 2004 of the optical light echo 
around V838 Mon using F435W as blue, F606W as green and F814W as red. Overlaid from left to right are the 24 (PSF subtracted - see text), 70 (PSF subtracted - see text) and 160 
fj,m contours. Contours are: 0.074 (1-cr), 0.11, 0.17, 0.26, 0.40, 0.61 mjy/arcsec 2 at 24 fj,m; 0.17 (3-<r), 0.28, 0.56, 1.23, 1.69, 2.26, 2.82, 3.38 mjy/arcsec 2 at 70 flia; and 0.52 (3-cr), 0.67, 
0.87, 1.13, 1.46, 1.89 mjy/arcsec 2 at 160 fim. The yellow line in the first image indicates a scale of 50". All images are oriented with north up and east to the left. 



